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CONTINUA 



P. VANi 



A. BEREZOVSKI* AND C. PAPENFUSS^ 



Abstract. Governing equations of dissipative generalized solid mechanics are 
derived by thermodynamic methods in the Piola-Kirchhoff frameworli using 
the Liu procedure. The isotropic small strain case is investigated in more 
detail. The connection to the Ginzburg-Landau type evolution, dual internal 
variables, and a thermodynamic generalization of the standard linear solid 
model of rheology is demonstrated. Specific examples are chosen to emphasize 
experimental confirmations and predictions beyond less general approaches. 



1. Introduction 

Conventional theories of continua do not provide the description of a microstruc- 
tural influence because material elements are considered as indistinct pieces of mat- 
ter. Generalized continuum theories (higher-order or higher-grade) are first exam- 
ples of what has been proposed to describe the macroscopic behavior of materials 
with inner material structure. Though their constitutive structure is restricted by 
the second law of thermodynamics, the contribution of small-scale events to entropy 
fluxes and sources is still not completely investigated. 

Governing equations of generalized continuum mechanics can be obtained by 
different ways. Three most widely accepted approaches are the following: the 
variational approach by Mindlin jl] , the microhomogenization procedure by Eringen 
an Suhubi [2 [3] , and the virtual power method by Germain [4] . 

Mindlin [1] derived governing equations in the small strain approximation with 
the help of a variational principle analogous to that in ideal elasticity. In his the- 
ory of microdcformations, Mindlin introduced kinetic and potential energy for both 
micro- an macro-displacements as well as a tensor characterizing micro-inertia. In 
the Mindlin theory, the potential energy density is a quadratic function of the 
macro-strain, the relative strain, and the micro-deformation gradient. With this 
variational foundation, the Mindlin theory is an idealized one, which does not in- 
clude any dissipation. 

In their approach to generalized continua, Eringen and Suhubi [21 |3] did not 
start from a variational principle. They obtained an evolution equation of the 
micro-strain extending mechanical concepts of inertia, stress, strain, and energy 
onto the microlevel and calculating velocity moments of the microevolution of the 
momentum. The zeroth moment of the mesoscopic momentum balance results in 
the macroscopic Cauchy equation and the first moment gives the evolution equation 
of the micromomentum. In their constitutive theory, Eringen and Suhubi consider 
the microdeformation and the gradient of the microdeformation as general internal 
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variables and calculate the entropy production accordingly. They obtained the 
following expression for the entropy production rA([2 (5.13)): 

TA ^Qidi logT + {Tij - fij)d.iVj + [ski - tu - aki)vki + 

Here qi is the heat flux, i.e. flux of the internal energy density, T is the temperature, 
Tij is the stress, p is the matter density, Sij is the relative stress, Uij is the microv- 
elocity gradient, /Xyj, is the double stress, the derivative of the free energy with 
respect to the microdeformation gradient. Tilde denotes the reversible, nondissipa- 
tive parts of the corresponding quantities. Further, (7^ = t^ — Sy is the reversible 
microstress, the derivative of the free energy with respect to the microdeformation. 
Equation ([T]) defines Sy, the reversible relative stress, too. The last term gives 
the entropy production due to additional internal variables ^a, a — 1, ■•"., and Ca 
denotes the related intensive variable, the derivative of the entropy density with 
respect to £,a ■ Here the notation of derivatives and indices follow Mindlin as far as 
it was possible. Remarkable is that the Eringen-Suhubi theory defines the entropy 
flux in its classical form Ji = qt/T. 

As we can see, the descriptions of generalized continua are weakly nonlocal from 
the beginning because the gradient of the microdeformation is introduced as a state 
variable. This is true also for the third approach that uses the principle of virtual 
power in order to derive the evolution equation of the microdeformation. The 
virtual power method is essentially dealing with statics. Dynamics in this method 
is introduced by an assumption that inertia is connected to virtual displacements 
directly [H [^ or with the help of dissipation potentials. In the first case, the 
dissipation requires separate considerations. In the second case, the dissipation is 
introduced together with inertial terms. 

Let us underline some common properties of above mentioned approaches: 

(1) The connection of new micro structural variables to mechanical effects has a 
kinematic background. However, the microdeformation can be originated in 
different structural changes (e.g. microcracking) , which is not necessarily 
connected to energy alterations due to the change of a Riemann geometry 
of the material manifold. 

(2) The evolution equations of micro structural variables are originated in me- 
chanics. Variational principles, moment series developement as well as 
virtual power with dissipation potentials are mechanical concepts. The 
dissipative effects are considered as secondary and frequently disregarded. 

(3) The entropy flux has the same form as that in the simplest case of Cauchy 
continua. 

(4) In all these approaches, it is customary to introduce internal variables in- 
dependently of the microdeformation, representing already identified struc- 
tural changes of continua, e.g. damage or cracking [6j[7]. 

Intuitively, it is natural to expect that the microstructure affects also dissipative 
phenomena. It is needed, therefore, to investigate in detail how dissipation effects 
are described by generalized continuum theories. 

The response of materials to external loads can be expressed explicitly as a 
functional [8]. On the other hand, additional internal variables can be introduced 
to define this functional in an implicit manner by means of their evolution equations. 
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It is well known, that variational principles may exist for a dissipative evolution, 
but they are not of the usual Hamiltonian kind. At best, they need to be modified 
because they do not work without any further ado [5J (TUl E]. Regarding the 
homogenization technics of Eringen and Suhubi for internal variables [2] , it should 
be noted that there is no primary microscopic candidate for an evolution equation, 
therefore their method cannot be applied beyond the kinematic determinations. 
Direct statistical or kinetic theory related calculations (as for example [12]) would 
require a particular microstructure and an interpretation of the internal variable. 
Finally, regarding the third treatment, the principle of virtual power is a mechanical 
concept, and internal variables (if any) are not related to any kind of spacial changes 
of continua. 

It is known, that there are two alternate approaches to obtain the evolution 
equations of internal variables directly. The thermodynamic approach is related to 
internal variables of state and the variational one is regarded to dynamic degrees 
of freedom, respectively |131 114] . In the thermodynamic approach, the entropy 
production is calculated considering additional constraints (e.g. balances), and the 
evolution of internal variables is determined as a part of the constitutive theory 
by means of the dissipation inequality. Evolution equations obtained in such a 
way contain typically only first order time derivatives |15] . The most important 
issues here are the stability of the weakly nonlocal extension (see critical remarks in 
[16|). the thermodynamic consistency, and seemingly missing boundary conditions. 
In the case of dynamic degrees of freedom, a Hamiltonian variational principle is 
applied to the nondissipative part of the evolution, and the dissipative contribution 
is calculated by dissipation potentials. 

Recently, it was shown that one can get a unified description of the two meth- 
ods introducing weakly nonlocal dual internal variables |17) . Then the exact and 
constructive exploitation of the entropy inequality, e.g. the Liu procedure, com- 
bined with an Onsagerian linear approximation of constitutive functions leads to 
completely solvable constraints. An essential ingredient of the approach of dual in- 
ternal variables is the observation that in the case of higher-order gradient theories, 
gradients of constraints of the entropy inequality are constraints on the constitutive 
state space [T8lp^l20]. 

It has been observed recently [JI] that the structure of the dual internal variable 
system of evolution equations in the nondissipative case corresponds exactly to the 
evolution equation in the Mindlin theory. In this case, one of the internal variables 
can be interpreted as the microdeformation and another as the conjugated momen- 
tum. It is worth to extend the dual internal variable approach to a broader class 
of materials taking into account dissipative effects. This suggests a more general 
procedure to construct evolution equations than in |17j . The corresponding pro- 
cedure is presented in the paper. It consists in the extension of the state space, 
the formulation of constraints, the application of the Liu procedure to the entropy 
inequality, the solution of obtained Liu equations, and the specification of a general 
form of evolution equations for internal variables following from linear relationships 
between thermodynamic fluxes and forces. The relation to linear viscoelasticity, to 
the pattern formation equations of the Ginzburg-Landau-type, and to the standard 
linear solid model is demonstrated in the small-strain approximation. The same ap- 
proximation is used to point out the extension of generalized continua descriptions 
onto dissipative materials and microstructural thermal effects. 
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2. Construction of evolution equations 

2.1. Balance laws. We start with the formulation of thermodynamic constraints 
for continuum mechanics with dual internal variables in the Piola-Kirchhoff frame- 
work (PK frame). This means that our treatment is only partially objective. 
The balances of momentum and energy can be represented as follows: 

goVi - djtij = 0, (2) 

goe + diQi ^ tijdjVi, (3) 

where po is the matter density, Vi is the velocity field, Uj is the first Piola-Kirchhoff 
stress, e is the specific internal energy and qi is the flux of internal energy density 
in the PK frame. The dot denotes the material time derivative, which is a partial 
time derivative on the material manifold, dj denotes the (material) space derivative. 
The Einstein summation rule is applied for repeated indices. Our index notation 
is abstract, does not refer to any particular system of coordinates, and denotes the 
tensorial degree and contractions in accordance with the traditional coordinate free 
treatment in continuum mechanics |22j . In the notation we use uniformly lowercase 
indices, i.e. do not distinguish between material and spatial indices and vectors and 
covectors. This way it is easier to follow calculations, and examples are considered 
in the small strain approximation where the differences are negligible. 

The balances are introduced without source terms, the momentum and the total 
energy are conserved, because source terms are irrelevant in a constitutive the- 
ory. The relation between the deformation gradient Fij and the velocity field is 
considered as a constraint: 

F,, - djv, = 0. (4) 

Evolution equations of the internal variables ipij and /3y are formally represented 
as 

Vij — Jiji Pij ~ 9iji (5) 

where the constitutive functions fij and gy depend on the whole set of state vari- 
ables. Here the notation of Mindlin was applied for the first internal variable, ipij, 
which is the microdeformation there [T] . 
The entropy inequality is given as follows: 

gos + a, J, > 0, (6) 

where s is the specific entropy and Ji is the fiux of the entropy density. 
It is assumed that constitutive functions 

'^ij : Qi J Jij : 9ij i ^i 'Ji-i K' ) 

are defined on the weakly nonlocal state space spanned by the following variables: 

djVi, Fij,dkF.ij,e, d.,e, 'il^ij,dk^ij,dki^ij,Pij,dkl3tj,dkiPij. (8) 

For the sake of simplicity, we consider a weakly nonlocal constitutive state space 
of the first order in the deformation gradient and in the internal energy, but of the 
second order in the internal variables, i.e., their second gradients are included. The 
velocity field is distinguished, because only its derivative is present in the constitu- 
tive state space. This assumption allows us to avoid velocity related problematic 
aspects of the material frame indifference (see e.g. ^5]). 
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2.2. Liu procedure. Balance of linear momentum ([2]), balance of internal energy 
(jSj, kinematic relation ([4]), and evolution equations of internal variables ([5]) are con- 
straints of the entropy inequality. Taking into account that the constitutive state 
space is weakly nonlocal of the second order in the internal variables, we should 
introduce additional constraints for derivatives of their evolution equations. We 
consider a derivative of a constraint as a new constraint in case of higher order 
weakly nonlocal constitutive state space. It is an important aspect for the devel- 
opment of correct thermodynamic conditions in weakly nonlocal thermodynamic 
theories [Ml [23: 



dk-fptj = dkfij dk$ij = dkgij- 



(9) 



Then we introduce Lagrange- Farkas multipliers A^, k, A^, Aij, Atjk, Bij, Bijk for con- 
straints ^, ©, (HI), ©I, 
imbalance is, therefore. 



L, ©2, &2, respectively. The constrained entropy 



QqS + djJj - Aj, (Pij - djvA - \i {goVi ~ dji^j) - k (poe + d,qi - t.jdjV,) + 



A 



J* 



B, 






' Jij 
9ij 



+ Akji ( dkipij + dkfij + 



(10) 



Bkji [dkPij - dkgi. 



>0. 



Liu equations are obtained as coefficients of higher derivatives after a straightfor- 
ward calculation: 
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ddjv.s = 
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eo dpi, s = Aij 
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dd,F,,s = 
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die 
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£•0 9^ij s = -A-ji 
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Qa do^^.^s = -Akji 
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£•0 dt3i^s = -Bji 
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Qo dd^PijS = -Bkji 
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(24) 
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dd,eJj - deS dg^e Qj " 
£'0 (^djpi^S dg-e gim 


QO ddjTpi^ 
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S do,e 
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(25) 
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dkim^ij : 9d,^-4,,jJk - deS ddi^^,,^ qk - Qo ddki>ap^ f^9imi/j.j fop- 

Qo do^p^^s da,^4,,^ gap = (26) 

dkimPij ■■ 5a,„^„- Jk - deS ddi„^[3,j qk ~ Qo da^p^pS da^^p^. fop- 
go da.fi^^s da,^0,, gop = (27) 

Liu equations ^, ^, ([15]), ^, ^, ^ and dH]) determine Lagrange-Farkas 
multipliers by corresponding entropy derivatives. The solution of Liu equations 
(HH), (dH), ((HI), Cni) and (1211) reduces the constitutive form of the entropy to the 
following one: s — s{Fij jejipijjdktpij, f3ij,dkf^ij)- The dissipation inequality then 
follows considering Liu equations (P5)) - (P7)) for the entropy flux. Together with 
the mentioned form of the entropy function, these equations can be solved and we 
obtain the entropy flux in the form: 

Jj = desq-t + goda,4,i„,s fim + Qoda^pi^s gim + J°- (28) 

Here the dependence of Jf is reduced similarly to that of the specific entropy 
Jf — Jf{Fij,e,tpij,dkipij,l3ij,dk(3ij). It should be noted that Eq. ([28]) is not the 
most general solution of corresponding Liu equations, because we did not consider 
symmetries of functions in Eqs. p^ . ([2^ and Q-^T^. For example, Eq. Q is 
obtained as multiplier of dkjVi, therefore only the symmetric part of Eq. Q must 
be zero, but we consider, however, solutions that are more restrictive. 
The dissipation inequality then follows as 

dkFij {dp^jJk - deS dp.j qk - Qa %>>,„« 9f,j fim - go da^i3i„,s dpi^ gim) + 
die {deJi - deS de qi - Qoda.^PimS de fim ~ go dd.Pi^s degim) + 

dk ^ij {d4,,jJk - deS d^,^. qk - goddk^im^ ^.^p^j fim " go da^/^i^s d^,^ gi^) + 
diki^ij {ddi^ijJk - deS (9a,^,j qk - go da^^p^^^s da^xp^^ fom — go dakPorr^^ (^dnf>,j 9om) + 

dkPij {dfi,.Jk - deS dfi^. qk - go da^Mn^s dp,^ fim - go da^Pi^s dp^. gi^) + 

dik Pi] {daiiSijJk - deS daip^^ qk ~ go ddtM>o„,s daip^^ fom ~ go da^Pa^^s daip^^ gom) + 

go dp^jS Fij +deS tijdiVj - go d^^.s fij ~ ga dp.^s g^ > 

(29) 

By substituting Eq. (j28p into the dissipation inequality we arrive at the following 
expression: 

dk (des) qk + deS U-jdiV-j + djVigodp.jS- 

- {9i,,jS- dkda^^ijs) go fij - {dp^^s - dkda^p.ys) gogtj + diJ^ > 0. 
Here deS ~ ^/T, and we may identify thermodynamic fluxes and forces as follows: 

5» f y j * + ^ {tjj + goTdp,^ s)^ d^vj - {d^,^s-dkda^4,,^s)^ eo/»j - 
- (9^,,s - dkda^p.ys) gogij + diJ° > 0, 



where t^j is the viscous stress. The entropy is a distinguished constitutive function, 
that fixes the static information of the system. All other constitutive functions are 
determined by means of entropy derivatives. 
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2.3. Evolution equations. The entropy production in tlie dissipation inequality 
is represented as a sum of products, and there is an undetermined constitutive 
function in each term multiplied by a given function of the constitutive state space. 
Therefore, it is straightforward to point out the simplest solution of the dissipation 
inequality assuming linear relationships between thermodynamic fluxes (terms with 
undetermined constitutive functions) and their multipliers, thermodynamic forces 
(see Table 1). 





Thermal 


Mechanical 


Internal 1 


Internal 2 


Fluxes 


Qi 


[Uj+Tgodp^^s] /T 


h 


ga 


Forces 


a.(i) 


diVj 


Qo {d^ijS—dkda^iPij s) 


Qo{dp.^s-dkda^p.^s) 



Table 1. Thermodynamic fluxes and forces 

The classical thermal interaction is vectorial, while other terms are tensorial. The 
mechanical term is responsible for viscoelasticity if no other terms are present. 
The last two terms constitutively determine evolution equations of second-order 
tensorial internal variables tpij and Pij. In isotropic materials, tensorial mechani- 
cal and internal variables may be coupled independently of the vectorial thermal 
constitutive function: 



goTdpi^s) 



Qt 


= A9,- 


2^ y 


= ^ {t^. 


A] 


= /.. 
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= 9 1] 





(32) 


ijkidkVi + Lij^iXki + L^j^.iYki 


(33) 


ijki^^kVi + Lijki^ki + LijkiYki 


(34) 


ijkidkVi + L^^ki^ki + L^jkiYki 


(35) 



So we have the three times three matrix of fourth order tensors if,,^; called conduc- 
tivity tensors. In isotropic materials they are characterized by three scalar material 
parameters each and can be represented for all /, J = 1, 2, 3 as follows: 

L"ki = li^S.kSji + li-'SaSjk + li-'5,j5ki. (36) 

where 5ij is the Kronecker delta. Therefore, there are 1+3x9 = 28 material conduc- 
tivity coefficients in isotropic media with dual internal variables. Representation 
(1361) can be decomposed into traceless symmetric, antisymmetric (deviatoric) and 
spherical parts, i.e., is equivalent to 

Lljki = s"S^kSji) + a^-'S,[kSji] + l"S,jSki, (37) 

where braces () denote the traceless symmetric part of the corresponding tensor in 
related indices ^^(^5^7) = {SikSji+SjkSii)/2—SijSki/3 and the rectangular parenthesis 
[ ] denotes the antisymmetric part as Si^k^ji] = {SikSji — Sjk5ii)/2. Therefore, s^'^ = 
li-' + li'^, a" = l{^ -li\ auAl" = {■ill^ + li-' + li'')/^. This kind of decomposition 
is instructive because symmetric, antisymmetric, and spherical second-order tensors 
are mutually orthogonal in the "double dot" product, i.e. taking the trace of their 
product. Therefore, constitutive equations ([55|) -([55 |) can be decomposed into three 
parts: five component traceless symmetric, three component antisymmetric, and 
one component spherical parts are independent. The spherical part is determined 
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as 

tl^ = l^^dkVk + l^^Xkk + ^"Ifefe, (38) 

i'kk = l^^dkVk + P^Xkk + l^^'Ykk, (39) 

Pkk = l^^dkVk + l^^Xkk + l^^Ykk- (40) 

In its turn, for the symmetric traceless, deviatoric part we have 

%) = S''^{^yl) + S''^(.,) + s'%^J) , (41) 

P(,^^ = s'^a^.z;,) + s32x<^,^ + s33y^,.^.^ (43) 

and the antisymmetric part is given in a tensorial form instead of the vectorial 
invariants: 

tm = a''d[,v^ + ai2X[,,] + fli^yj^^.j ^ (44) 

ip[,j^ = a^^d[,v^ + a22X[y] + a'^^Y[,j] , (45) 

%] = a^^di.vi] + a32X[,,] + a^^yj^^.j . (45) 

Therefore, conductivity matrices of spherical, traceless symmetric, and antisym- 
metric components are 

lU ^ pi ;22 ;23 \ /J ^ ^21 ^22 ^23 ^ „/,/ , 

\/31 ;32 ;33y \^^31 ^32 ^33^ 

(47) 

The second law requires that symmetric parts of 3x3 conductivity matrices l^'^ , 
s^'^ , and a^"^ are positive definite. 

With this requirement the thermodynamically consistent construction of con- 
tinuum mechanical theory with second order tensorial dual internal variables 
is complete. Balances of momentum and energy ^, ^ and evolution equa- 
tions of the internal variables ([5]) together with isotropic constitutive functions 
(|32l) - (j35p form a complete system, if a particular form of the entropy function 
s {F^j , e, V-jj , dkipij , /3jj , dkPij ) is given. 

2.4. Boundary conditions. It is remarkable that one can get natural boundary 
conditions for internal variables with the requirement that the intrinsic part of 
entropy flux (j28l) is zero at the boundary. There are three basic possibilities: 

(1) Combined condition. In this case dg.^^^s and dg-pi^s are orthogonal to fim 
and gim- 

(2) No change condition. Internal variables do not change at the boundary, ipij 
and I3ij are zero. 

(3) Gradient condition, tidg.^^^s = Oim and tidg-^^^^s — Oim for any vector ti 
that is parallel to the boundary. In the case of a quadratic dependence on 
gradients, the condition is that gradients of internal variables are orthogonal 
to the boundary. 



THERMODYNAMIC APPROACH TO GENERALIZED CONTINUA 9 

2.5. Remarks on reciprocity relations. Classical irreversible thermodynamics 
requires special conductivity matrices obtained by reciprocity relations. In the case 
of symmetric matrices, the reciprocity relations are of Onsager type |26[|27) . if the 
matrices are antisymmetric, then they are of Casimir type [28 . These restrictions 
are justified by arguments from statistical physics requiring a microscopic inter- 
pretation of thermodynamic variables. If these variables are even or odd functions 
of the microscopic velocities, then the conductivity matrix is symmetric or anti- 
symmetric, respectively. However, in our case one cannot specify the conductivity 
matrices, the conditions of Onsager or Casimir type reciprocity relations are not 
fulfilled. In the following we will show some possible interpretations of the dual 
internal variables, and we will see that the most straightforward interpretation, 
namely, the micromorphic generalized mechanics when the internal variables are of 
deformation and deformation rate type, requires general forms of the conductivity 
matrices. Moreover, the background of internal variables in continuum mechan- 
ics is originated in structural changes in a material, but not on the microscopic, 
atomic, or mesoscopic level [17]. Therefore we simply do not see any reasons to 
further specialize the theory and we keep our treatment universal, independent of 
microscopic or mesoscopic interpretations. 

Idealized theories are characterized by a nondissipative behavior. In the case of 
heat conduction, the heat conduction coefficient A is zero and, therefore, the heat 
flux, the flux of the internal energy density, is zero, too. In the case of pure mechan- 
ical interaction without internal variables, the idealized theory is characterized by 
zero viscosities, the bulk viscosity l^^ ^ r/y = and the shear viscosity s^^ = 2ri = Q 
are zero, together with the condition a^^ = 0. Therefore, the viscous stress is zero, 
the stress in the momentum balance is the static stress t*-' = —gTdpiS. However, 
antisymmetric terms of conductivity matrices do not produce entropy, therefore in 
our case the coupling between different terms may result in nondissipative trans- 
port, too. 

3. Spatial representation 

In this section, we represent governing equations in spatial framework and at the 
same time introduce the small strain approximation. Technical details are given in 
the Appendix. 

We can transform balances of internal energy and momentum ([2]) , ^ into spatial 
form due to the Nanson theorem, to obtain 

ge + diqi^ UjdjVi, (48) 

gv, ~ djUj = 0. (49) 

Here g ~ go is the density in the actual configuration, tij is the Cauchy stress 
tensor, and Qi is the spatial fiux of the internal energy, the upper dot denotes the 
substantial time derivative, e = dtc + v^d^e. Spatial and material forms of these 
balances are identical in the small strain approximation. 

The spatial form of evolution equations of internal variables is the following: 



The symbol ^ denotes the spatial form of the material time derivative of a second- 
order tensor. The material time derivative is an upper convected one, because 
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internal variables are defined on the material manifold, supposedly independent of 
the motion of the continuum (see Appendix): 







ipi 



'ij- -vtj - dkVitpk] - dkVjiJjjk, (51) 

i.e., spatial and material forms of constitutive functions for internal variables are 
the same in the small strain approximation if velocity gradients are negligible. The 
spatial form of the entropy flux is represented as 

Ji = deSQi + gda.^^^s fim + Qdo.Pi^s gim + J^ ■ (52) 

We do not deal here with any interpretation of the extra entropy flux j':^ and, 
therefore, we assume that it is equal to zero. This requirement is not necessary, 
and there are phenomena which can be modelled only with the help of nonzero 
extra entropy flux [3D1 [211 [32] . 

Finally, the spatial form of the entropy production in the small strain approxi- 
mation is the following: 



dA^ 



qi 






gTdeij s) k 



(53) 



- {d^,,s - dkda^4,,js) gfij ~ {dp^.s - dkdo^p.^s) ggij > 0. 

Here e^j is the Cauchy strain. 

The original dependency of the entropy function on the deformation gradient 
Fij turns to a dependency on e^ because of the small strain approximation. More- 
over, it is convenient to introduce the specific free energy function instead of the 
specific entropy as a thermodynamic potential. In our case, the expression of the 
corresponding partial Legendre transformation is: 

w{T, ey , ?Aij, c'fcV'y , ftj, 9fcftj) = e - Ts{e, eij,^pij,dkipij, I3ij,dkl3ij), (54) 

Therefore, 

ds 1 dw 

-— = — and — — 
de T dT 

keeping fixed all other variables. Other partial derivatives of w and s are related 
as follows: 

ds 
de,. 



-s, 



(55) 
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1 dw 

T ddktp^j 
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ddkPij 
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1 dw 



T 5/3y 
1 dw 



TddkPr 



(56) 



All not indicated variables are kept fixed in the partial derivatives. In terms of the 
free energy function, we can transform entropy production (j53p to the form 

Ta = di InT [qi + gfijdg^^^.w -\- ggijdg^p.^w) + [Uj - gd^^.w) e^H- 

(57) 
+ [^i^^JW - dk{do^i,,jW)) gf^j + [dp^^w - dk{dd^j3,^w)) gg^j > 0. 

It is worth to introduce thermodynamic forces and fluxes in the small strain ap- 
proximation according to Eq. (|57l) : 





Thermal 


Mechanical 


Internal 1 


Internal 2 


Fluxes 


9fc = qk + 




Qfij 


Qijij 


Forces 


dklnT 


iij 




dfi-- w - dk dg^ 11. . w 



= Xdk\nT, 


(58) 


L]]ki^ki + L]]kiXki + Llfki^ku 


(59) 


Lfjki^ki + LlfkiXki + i4%iYkU 


(60) 


L%ihi + L^fki^ki + Lf^kiYki- 


(61) 
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Table 2. Thermodynamic fluxes and forces in the small strain approximation 
It is easy to see that thermodynamic forces and fluxes in the small strain approxi- 
mation are very similar to those in the PK frame. The most important difference 
is the regrouping of terms which are proportional to the temperature gradient. 
This representation of forces and fluxes is more convenient for the separation (or 
coupling) of thermodynamic and mechanical parts of the entropy production, espe- 
cially in the case of thermal stresses. The solution of dissipation inequality ([57)) is 
provided by conductivity equations: 

qk = qk + 9fijddk4>,jW + ggijdd^p.^w 

iij = Uj - gde^.w 


Q -fPrj = Qfij 



Q Ptj = ggij 

In the case of isotropic materials, a decomposition of conductivity matrices can 
be performed introducing symmetric, antisymmetric and spherical parts of corre- 
sponding tensorial forces and fluxes similarly to previous case: 

U/ki = s"S^kSji) + a"S,[kSji] + l"S^jSki. (62) 

In the following, we treat different special cases, all of them with constant conduc- 
tivity matrices. 

4. Examples 

4.1. Linear viscoelasticity, relaxation, and Ginzburg-Landau equation. 

We start with the simplest case when there is not any coupling between evolution 
equations (|58 p - (|61|) . In this case, the balance equation of the linear momentum 
and evolution equations for internal variables are independent. This means that 
the conductivity hypermatrix is diagonal L^^ — L^^ = L^^ — L^^ — L'^'^ — L^^ — 0, 
and we assume that the heat conduction coefficient A is zero. The free energy is ad- 
ditively decomposed into parts which are dependent on ipij , /3ij , and eij separately. 
Dissipation inequality (fST)) reduces to 

Tcr = {Uj - Qde^^w) hj + 

(o-j) 

+ [d-4>,jW - dk{dd^i,i^w)) Qf,j + [dfi^.w - dkidg^Pi^w)) gg,j > 0, 

and we see that terms related to internal variables are completely similar. This 
means that in the absence of couplings it is enough to analyze only one of them. 
For the viscous part of the stress we have, therefore, 

Uj - tf; = t^, = l^^hkkj + s"e(,,) . (64) 

where the elastic stress is introduced as usually 

g-^ — =tf°' = g\ekk5^J+ g^iitij, (65) 

l^^ — Vv corresponds to the bulk viscosity and s^^ — r/ ia the shear viscosity. 
All other coefficients in conductivity matrices are zero due to the absence of any 
coupling. The evolution equation for the internal variable ipij results in: 

.„„ 

e i^tj = L,ji^ {di,i^w - dkdd^^.^w) . (66) 
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Assuming the isotropy of the conductivity tensor L^^ , we can decompose Eq. (|66p 
into six independent evolution equations for the spherical, symmetric traceless 
and antisymmetric parts of the internal variable tensor. These evolution equa- 
tions for the internal variable tpij give the generalization of the Ginzburg-Landau- 
Khalatnikov equation, introduced first as a scalar equation in the case of supercon- 
ductors [34l|35]. For each free energy functional, the structure of such equations is 
universal and widely used with different thermodynamic arguments [36^ [37^ i38| 139) . 
The Ginzburg-Landau equation was derived by pure thermodynamic arguments 
as the evolution equation for a second-order weakly nonlocal internal variable in 
[24] (see also [21] )• That derivation shows a universal character of the Ginzburg- 
Landau equation: the second law requires an evolution equation of this form for 
an internal variable without any other constraints independently of the microscopic 
background. 

The symmetric traceless form part of Eq. (j66p gives the de Gennes-Landau 
theory of liquid crystals, if a suitable quadratic dependence of the free energy on 
'ip{ij) is introduced |40j . This very particular example of the proposed approach 
shows the richness of the mathematical structure and its physical interpretation. 

4.2. Generalized standard linear solid. In this section we consider a coupling 
between a single tensorial internal variable and deformation, which results in the 
thermodynamic theory of rheology [411 HI]. Therefore, we still assume zero heat 
conduction coefficient and L^^ — L^^ = L^^ = L^^ — 0. However, L^"^ ^ and 
LP'^ 7^ anymore. Let us further reduce the treatment by introducing a local theory 
for the single internal variable with the simplest quadratic dependence of the free 
energy on the internal variable without its gradients. In this case, the free energy 
can be written as 



w(ey, V'y) = :t4 + /^e-ye«j + 



(67) 



"^V'ii + ~^'4'ij''Pij + "^ViJ-A-Ji -r yi^iz-^jj r yzyvij "r f/jjycij. 



Thermodynamic stability requires the convexity of the free energy, hence the in- 
equalities follow: 

3Ai + 2A > 0, Ai > 0, ^2 > 0, 62 + 63 > 0, 3&1 + 2&2 > 0, 
M62 + 63) - 5? > 0, (3A + 2^)(36i + fe2 + fo3)-(3gi+2g2)'>0. ^^^^ 

The entropy production contains two tensorial terms that may be coupled: 

Ta = {Uj - gd.^^w) hj - d^^^wgP > 0, (69) 

In isotropic materials, spherical, deviatoric and antisymmetric parts of the ten- 
sors are independent. For the deviatoric, symmetric traceless part we obtain the 
following conductivity equations from Eqs. (p7)) and (|69p : 

t{zj) - 2£i(/xe(,j) + g2^{ij)) = s'^^iij) - s^^ ((62 + &3)^(y) + '^92<^{ij)) (70) 

S V'(.y) = s'^^^^■J) - s^ ((&2 + h-i)iPt^,j) + 2(72e(»j)) • (71) 

Equation ([70]) is a constitutive equation for the deviatoric part of the stress i^^j) , 
and Eq. ([7T|) is the evolution equation of '4'{ij) ■ Moreover, the second law requires 
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that the symmetric part of the conductivity matrix is positive definite, therefore, 

'12 _ „21N 2 
„1 \ n „2 \ A „1„2 „12„21 ' * * 



s^>0, s^>0, s's'~s''s''~ r ^ I >0. (72) 

The role of the internal variable may be better understood if we eliminate it from 

Eqs. ([70|) ~([7T 1) . Taking the material time derivative of Eq. (fTO]) and substituting 



i^Hj) from Eq. ([7T|) and V'(ij) from Eq. (|70p into the obtained form, the following 

relation follows: 



s'e/,.^ + p [{b2 + &3)(2m + (s^s' - s''s'') + (s'' - si')2g2)] 



(y) + F L'^'^2 + osn^M + l« * -* * ^ + 1* -* J^y2jJ e(,j)- 

2„2^^7 , , N,, o„2n 



2sV((&2+&3)A*-V)efe) (73) 

The positive definiteness of the free energy w requires that the coefficient (62 + 
&3)/i — 2g^ is non- negative and Eq. ([75|) can be transformed to 

where 

ps^{b2 + b3) 
y = P [(62 + b3){2p + (s\s2 - fii^s^i) + (,s2i - ,si2)2g2)] T. (75) 



Constitutive relation (|74|) reduces to the Zener body of rheology or standard linear 
solid [43l |44] , if t^ = and the nonlinear part of the time derivative is neglected. 
This rheological model is widely used in different fields from biology [45] to en- 
gineering |46) . The complete form of this constitutive relation (with r^ 7^ 0) is 
called the inertial Poynting- Thomson body. The proposed thermodynamic model 
has remarkable properties which are distinctive in comparison with more intuitive 
approaches: 

(1) Neither the inertial term beyond the standard model, nor the coupled 
volumetric-deviatoric effect is neglectable in general, both are important, 
e.g. in experimental rock rheology [1711351112]! where an inertial Poynting- 
Thomson body is used for both deviatoric and spherical parts of the defor- 
mation in order to get a good agreement with experimental data. 

(2) Reciprocity has not been required in the proposed approach. The same 
is true in the simplest case of a standard linear solid body, otherwise the 
second law of thermodynamics contradicts to observations [41] . 

(3) The complete form of constitutive equation ([7n |) - (|7ip is preferable instead of 
((7^ if dynamical problems with the coupled balance of momentum ^ are 
needed to be solved. This is an important advantage of the thermodynamic 
approach. 

4.3. Dual internal variables. Now we consider dual (coupled) tensorial internal 
variables which are independent of mechanical and thermal interactions (A = 0, 
L^^ = L^^ = L^^ = L^^ = L^^ — 0). The evolution equations of the coupled 
tensorial internal variables again follow from the dissipation inequality. In the 
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small strain approximation and with small velocity gradients, Eqs. (|60p - (l6ip are 
simplified to 

i'lj = ^ijki^ki + L^j^iYu, (76) 

Pij = Lij^iXki + L.^jkiYki. (77) 

Here Xij = d^.^w ~ dk (ddktpijw), and Yij = dp^^w - dk {ddkPijW). Equations ([76])- 
([77]) are independent of balances of linear momentum and energy, if the free energy 
can be decomposed into a sum of functions depending on two sets of variables, e, Cij 
and tpijjPij separately, and if the objective time derivatives can be substituted by 
the substantial time derivative, that is the nonlinear terms are neglectable in the 
upper convected derivatives. 

In the case of isotropic materials, tensorial equations ([7S)) -([77 | are decomposed 
into a spherical, deviatoric and antisymmetric parts with scalar coefhcients. The 
remarkable difference between the evolution described by a Ginzburg-Landau-type 
equation based on a single internal variable and this dual structure becomes appar- 
ent after the separation of the symmetric and antisymmetric parts of the particular 
decompositions. 

As an example, we consider completely decoupled deviatoric evolution equations. 
Then conductivity matrices are two-dimensional and evolution equations are the 
following: 

'4>{ij) = si^(ii) + ('"^ + °)Y{ij) ' C^s) 

P{ij) = {s - a)X(^ij) + S2>'(ij) ■ (79) 

The above evolution equations are decomposed into a symmetric part, which repre- 
sents a dissipative evolution and therefore produces entropy, and the antisymmetric 
part that does not produce entropy and represents a non-dissipative part of the evo- 
lution. The role of the non-dissipative part can be better understood with the help 
of the following free energy function, where we assume a local theory for Puj) '■ 

Q 

W{f3{ij),1p{ij),dkll'{ij)) = ■^P{t3)l3(tj) +Wh{lp(^^j)) +Wg{^ki){^j)) = 

■^ 7 /» r p 

i^P{i3)P(i]) + 2^(y>^{y> + y^*;V'fe>5feV'(jj) + -i^dki^(ij)diil^(jk) + -jdkip{tk)djip{ij). 

(80) 

Here Wh is the homogeneous, local part and Wg is the gradient dependent, weakly 
nonlocal part of the free energy related to the variable ip{ij) ■ The second line of Eq. 
(|80p shows a particular quadratic form of corresponding functions. If /2 = /s = 
0, then we obtain the usual second gradient theory. The non-dissipative part of 
evolution equations (|78 |) -([79 |) has the form 



V'(y) = ac/3(^,j), (81) 

P{^3) = -a (9^(,,>w^ - dk{do,i;^,^^Wg)) (82) 

It is easy to eliminate /3(iA:) and obtain an evolution equation for 'ip(^ik) , which is 
second order in time: 

-2-^^m + d^^.k)Ws ~ dk (9a,v<,fc)W^s) = (83) 
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The last equation can be considered as the Euler-Lagrange equation of the La- 
grangian 

L{-4^{ik),i'{ik)) = ^^'/'{«fc>'0{ifc> - w(V'(ife),9fcV{»fc>)- (84) 

It is remarkable that natural thermodynamic boundary conditions of the zero en- 
tropy flux requirement correspond exactly to natural boundary conditions of the 
variational principle. 

4.4. Dissipative generalized continua. It was already observed by Berezovski, 
Engelbrecht and Maugin [2T| that generalized thermomechanics of solids is a partic- 
ular case of the dual internal variables theory. They observed that evolution equa- 
tions of the non-dissipative theory of Mindlin correspond exactly to non-dissipative 
evolution equations of the dual internal variables theory. However, they did not 
perform a complete thermodynamic analysis and, therefore, their observation is 
restricted to the idealized non-dissipative case. 

The corresponding thermodynamic analysis was performed in the second section 
of the present paper exploiting the dissipation inequality. The coupled evolution 
equations are represented in the form of linear conductivity equations. If one of the 
internal variables is interpreted as a microdeformation, then our calculations are 
to be considered as a pure thermodynamic derivation of a generalized dissipative 
continuum theory. The particular example of the theory of Mindlin arises then 
under the following conditions: 

(1) No thermal and viscous dissipation. 

(2) Pure antisymmetric coupling between internal variables. 

(3) Quadratic free energy function. 

Another generalized continuum theory has been introduced by Eringen and 
Suhubi [2 . To compare the entropy production, let us consider for simplicity a 
continuum without additional internal variables in the small strain approximation 
and uniform temperature field. Then the entropy production in the Eringen-Suhubi 
theory ([T]) lack the first and the last terms and can be written with our notation 
as 

Ta = {Uj - gde^jw) e^-, + {sij - nj - gd^,-w) ^ij + {^Hjk - dd^-^i.^w) dki^ij > 0. 

(85) 

Here the internal variable ipij is identified with the microdeformation gradient x'ij 
of the Eringen-Suhubi theory and therefore the material time derivative is the sub- 
stantial time derivative due to the deformation interpretation. For the comparison, 
let us repeat here our entropy production ([63)) : 

(86) 







The difference in the entropy flux in the two theories and the different concept of 
constitutive quantities determine the diversity in the entropy production. We have 
introduced the evolution equations of internal variables as constitutive relations to 
be determined from the entropy inequality. Eringen and Suhubi 2 simply indicated 
the form of the dissipation inequality following from their definition of stresses. 
Regarding the exploitation of the second law one should observe the following: 
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- The microniomcntum balance and the evolution equation of tpij are not 
constructed from the dissipation inequality. 

- The inicromomentuni balance is not used as a constraint for the entropy 
inequality in the Eringen-Suhubi derivation. However, it is implicitly con- 
sidered during the application of the Coleman- Noll procedure assuming that 
the multipliers of the time derivatives should be zero. 

- The entropy flux is not an arbitrary constitutive function, but it is restricted 
to the classical Ji — Qi/T in the case of Eringen and Suhubi. 

- The constitutive state space is not weakly nonlocal and not fixed in the 
Eringen-Suhubi derivation. 

As we have shown above, not only the basic assumptions, but also the final 
equations of the Eringen-Suhubi theory are particular and can be obtained from 
our generalized approach if several dissipative terms are neglected. Virtual power 
approaches also introduce an entropy production that have a similar form and 
similar limitations as the Eringen-Suhubi approach has (see e.g. [50ll51) ). 



4.5. Heat conduction and weakly nonlocal internal variables: Microtem- 
perature. Finally, let us consider the case with a non-zero heat conduction coef- 
ficient A. Neglecting the viscosity infiuence (L" = L^^ = L" = L^^ = L^^ = 0), 
we chose the free energy dependence on internal variables in the form of Eq. (I5n|) , 
but with the reduced conductivity matrix for deviatoric evolution equations 

V'(y) = aY^.j) , (87) 



which corresponds to the choice s = si = in Eqs. ()78|) -([79 |) . The free energy 
dependence (|80|) allows to represent the evolution equations in the form 

'0(y) = ac/3^,j) , (89) 

/3(jj) = -a (d^^..^Ws - dk{dd^^^..^Wg)) + S2cl3^,j), (90) 

which can be reduced to the single second-order evolution equation for the primary 
internal variable '0(y) 

-^^^m - -f ^('fc) + ^i'i^u)Ws - dk {^9^^^,^^wg) = 0, (91) 

which is similar to the Jeffreys type modification of the Maxwell- Cattaneo-Vernotte 
equation [55]. 

The thermal part of the dissipation inequality is satisfied by the modified Fourier 
law that follows from Eq. ((58|) 



qk + QhjddkiPijW -t- ggijddki3,jW = Xdk InT. (92) 

As it was shown [21] on the example of one-dimensional thermoelasticity, the pri- 
mary internal variable tp(ij} can be interpreted in this case as a microtemperature. 
In this context, it is understood as a fluctuation of the macrotemperature due to the 
influence of the existing microstructure. The solution of the equations shows that 
influence of microtemperature may result in a wavelike propagation of temperature 
if the corresponding damping effects are small [21] . 
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5. Summary and discussion 

The paper is devoted to the answer of the fohowing question: How could we 
obtain evolution equations of physical quantities, about which we do not know any- 
thing, i.e., only general principles can be considered? There are essentially two 
basic approaches. The first one postulates a variational principle of Hamiltonian 
type coming from mechanics. In this case dissipation is something additional to 
the non-dissipative basic mechanical evolution. The second approach is coming 
from thermodynamics: one can assume that the evolution of new variables is not 
exception from the second law and generate their governing differential equations 
accordingly. This is the situation in the case of internal variables in general as it 
was summarized by Maugin and Muschik [13[ 114] . The two approaches can be gen- 
eralized. Thermodynamic principles and dual internal variables in the framework 
of a second order weakly nonlocal theory give a straightforward and simple way of 
the generalization [17^ . 

As a result, the thermodynamic consistency of continuum mechanics with dual 
tensorial internal variables was analyzed in the present paper by the Liu approach 
to the exploitation of the second law in the Piola-Kirchhoff framework. Then local 
evolution was considered in isotropic materials in the small strain approximation. 
The entropy production was calculated and thermodynamic forces and fluxes were 
identified. Then a quadratic free energy and linear conductivity relations closed 
the system of equations. The final evolution equations in a non-dissipative case 
are equivalent of those for micromorphic continua, therefore, the thermodynamic 
method gave a dissipative extension of the original Mindlin theory. 

We have given several particular examples that arise as special cases of the gen- 
eral theory. Our goal was only partially a justification, but also the identification of 
the most important differences from other theories and the interpretation of some 
qualitative predictions of our approach. We have seen that phenomena of microtem- 
perature, sophisticated couplings in generalized rheology, and special properties of 
the dissipative extension without assuming reciprocity relations, all are open for 
experimental testing. We think that this approach is essential in the case of gener- 
alized continua, where additional coefficients are considered hardly measurable. In 
this respect we have analyzed the limitations of the dissipation in other classical 
approaches, in particular in the Eringen-Suhubi theory, which is one of the most 
developed from this point of view. We have seen, that due to the restrictive start- 
ing assumptions (mechanical interpretation of the internal variable, locality, special 
entropy flux, etc.) the considered dissipation is extremely limited. For example, 
under the traditional approaches one cannot recover neither the Ginzburg-Landau 
equation, nor simple viscoelasticity. 

It is remarkable, that the finite deformation part of our approach shares short- 
comings of the Piola-Kirchhoff framework. There are indications that the require- 
ment of objectivity and material frame indifference are not treated properly in 
this case [531 154] . Beyond the reservations of using material manifolds in general, 
it is also remarkable that the exclusion of velocity field Vi from the constitutive 
state space is not necessary [23] 155] . We expect several interesting phenomena by 
the analysis of higher-order nonlocality at the mechanical and thermal side. Here 
the comparison to phase field approaches looks like a promising direction (see e.g. 

[SiEZlEH]). 
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6. Appendix 

In this section we shortly derive the material time derivatives, introduce the 
small strain approximation and describe the transformation of the balances between 
a Piola-Kirchhoff and local frameworks. Here we distinguish between contra and 
covariant as well as between material manifold and space-time vectors and tensor 
components. The covariant and contravariant vectors are denoted by lower and 
upper indices, the space vector and tensor components are denoted by minuscules, 
and the vector and tensor components at the material manifold by capital letters. 
We assume here that the reference configuration is relaxed, stress free, therefore 
the transformation between the material and spatial descriptions is standard (see 
e.g. [33J- For a more detailed kinematics, considering general bodies, see [53] . 

The material vectors are denoted by X', the spatial ones by x*. Therefore, the 
deformation gradient, the material manifold derivative of the motion, x*(i,X'^), is 
given as: 

F\, = djx'. (93) 

The transformation between material and spatial vectors and covectors is the fol- 
lowing: 

a-' = {F-yy, a' = F\ja'', (94) 

bj = F'jh, h = {F-y,b.j (95) 

In particular, the transformation of space derivatives follows the lower indexed 
covector rule: 

dj = FV5„ 9,; = {F-y,dj. (96) 

The summation over repeated indices is still assumed. 

6.1. Material time derivatives. The transformation of time derivatives is dif- 
ferent for quantities with different tensorial character. For scalars, the partial time 
derivative on the material manifold, 9^, is the substantial derivative for local quan- 
tities. We use the convenient dot notation for local quantities that corresponds to 
the partial time derivative on the material manifold for scalars, for the velocity and 
for acceleration fields: 

dta{t, X') = a = dta{t, x') + v'^^a{t, x'), (97) 

dtx'^v\ dux'=dtv'^v\ (98) 

The material time derivative of internal variables with various tensorial character 
differ from each other. Here we give the calculation for the second-order tensorial 
variable, '0'"' • According to the definition of material time derivatives, the par- 
tial time derivative on the material manifold expressed by spatial fields gives the 
material time derivative: 

dt^" =dt {{F-^)\{F-y^^^') 

={F-ydF-y,i^'' - {F-y^dtF\{F~')\{F-y^^^^ 
-{F-y,{F-y^dtF\{F-^)^^r 

= {F-y,{F-% {^'^ - dkv'^"' - duv^i^'') . (99) 
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Therefore, the spatial form of the abovementioned formula, the material time de- 
rivative of the tensor, is given as 


^^J^ F^^pJjQ^^iJ ^ ^y _ dkV'^^^ - dkV^^'^. (100) 

Here we used the kinematic relation for the spatial velocity gradient and the time 
derivative of the deformation gradient djv"^ — dtF''j{F~^y',. In case of cotensors or 
mixed tensors, the spatial form of the material time derivative is different. 



6.2. Small strains. The small spatial strains are defined with the left Cauchy- 
Green deformation A^^ — F'^jF/ , as 

e'3 :=^{A'^ -S'^). (101) 

This choice is the best considering the requirement of objectivity [53]. The material 
time derivative of the strain in the small strain approximation is the symmetric part 
of the velocity gradient 

It should be noted that in the small strain approximation one may obtain identical 
results starting from different deformation concepts. 

The spatial form of the second material time derivative of the strain is particular 

e^J' = ii^J « 9('{,') + dkv'dkv^. (103) 

Here we have used that div'A'-i = F}{F-^)'\F''^fJ . 

6.3. Spatial balances. The relation between the local density g and the material 
density qq is 

= go 
^ detF' 
where det F is the determinant of F*^ 

The local and Piola-Kirchhoff forms of the heat flux and the stress are, respec- 
tively, 

q'^ {det Fy'^F'jq'', q' = dciF{F-^)/ q\ (104) 

f^ = {deiF)-^F\t'^\ t'^ = det F{F~y,,t''^. (105) 

For the transformation of basic balances, the Nanson theorem is essential. It can 
be written as: 

dj{detF{F-^y{)=Q,. (106) 

The proof is straightforward, when considering that the derivative of the determi- 
nant is 9,/(detF) = det F(F~^)-^(9ji^'^ and the derivative of the inverse deforma- 
tion gradient is dj{F-^y^ = -{F-^)\djF^j^{F-^)^j. 

Then the transformation of the balance of internal energy follows by substituting 
the definitions: 

Qoe + dKq^ = det Fge 4- det F{F^^Y.diq^ = 

, I- (107) 

= deti^(6ie + djq^) = t^^djiP = det F t^^ dkV^ . 

Similarly, the balance of linear momentum can be easily obtained: 

Q^v' + dKt'^' = ... = det F{ge + dkt''') = 0*. (108) 
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Therefore, we obtain usual local balances of internal energy and momentum (|48p 
and (I49p without approximations. 
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